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I. INTRODUCTION 


The determination of the frequency-dependent permeability tensor of a mag¬ 
netic unsaturated gyrotropic medium such as ferrites materials is a long stand¬ 
ing problerrti^-. It has been demonstrated useful to describe precisely this tensor 
in frequency to obtain novel radiation and scattering characteristics of microstrip 
antennas^. To the physical point of view in dense polycrystalline materials, the 
frequency behavior of the permeability in the range from lOOkHz to 20GHz is under¬ 
stood as a superposition of magnetization changes due to domain-walls motion and 
coherent rotations of magnetic domainsi*^. However, Grimes^*^ has suggested that 
the observed spectra may be also interpreted as multiple scattering of electromag¬ 
netic waves in these polycrystalline materials by both random magnetic and dielectric 
homogeneous spheres. Unfortunately this deduced interpretation has not been tested 
when magnetic materials, subject to external dc uniform magnetic held present an 
hysteresis behavior on the magnetization and changes the derived permeability. A 
convergence of these two alternate mechanisms was hrst recognized by Schlomann^ 
who’s derived a relationship between the isotropic permeability in a completely de¬ 
magnetized state by computing the scattering of magnetic helds in the static limit 
by out-of-phase and concentrical cylinders of gyromagnetic single-domains. The ex¬ 
tension of this idea to unsaturated magnetic states was conducted by Bouchaud 
and Zerah4 but the bridge between the local remanent magnetization, included as 
fractions of magnetic volumes of scatterers and the observed average distribution of 
magnetization directions of domains in the tridimensional space reminded elusive. 
To get into account the anisotropic nature of this scattering problem, Stroud^ was 
pioneer to generalize an old effective-medium approximation (EMA) for the con¬ 
ductivity tensor of a randomly inhomogeneous medium to treat materials consisting 
of crystallites of arbitrary shape and conductivity tensors of arbitrary symmetry. 
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Derivation of non diagonal permeability tensors for magnetized granular composites 
has been then reportedly and several models for the calculation of complex perme¬ 
ability of magnetic composite materials have been proposedl^y^. By using the EMA, 
even the influence of porosity induced by non-magnetic inclusions in magnetized het¬ 
erogeneous materials has been computediy. However, even in the quasi-static limit 
of frequency, the problem of taking into account dynamically the magnetic multiple¬ 
scattering of both do main-walls and single domains simultaneously in a consistent, 
anisotropic formulation is not achieved yet. This is the main purpose of this paper. 

II. EFFECTIVE MEDIUM APPROXIMATION FOR 
UNSATURATED FERRITES 

Because electrical non-conducting media are considered, the ferrites are treated 
as uniform from the dielectric point of view. The amplitude of an ac magnetic 
held applied is sufficiently low to let hrst the transverse permeability tensor being 
independent of it and to generate a small local transverse oscillating magnetic held 
only. 

Saturated ferrites exhibit anisotropic permeability described by the well known 
Polder tensor—, resulting in the nonreciprocal behavior of microwaves into magne¬ 
tized ferrite used for the design of circulators and isolatorsl^. As any reference axis 
may be chosen to project the third component of an arbitrary magnetization vec¬ 
tor, the ; 2 -axis is then taken for that. So in a cartesian frame and for any arbitrary 
direction, the magnetization vector is represented by only 2 angles, a polar one 6 
and an azimutal one (j). Then in this arbitrary frame, the gyromagnetic permeability 
tensor can be expressed in a general forn>iyy^. Such saturated ferrite is thus char¬ 
acterized by only three circular permeabilities which come from the response of the 
magnetization to a rotating ac magnetic held. This corresponds to diagonalize the 


3 


general form of the gyromagnetic permeability tensor as first recognized by Tyras 
a long time ago^^. It is always possible to select a frame that puts the saturation 
magnetization vector along the ^-axis and in that case, one of the third component 
of the permeability is equal to 1 and the ferrite is characterized by only 2 remaining 
permeabilities. 

For an unsaturated polycristalline ferrite, the spectral value of the 3 eigenvalues 
as a function of a given reduced remanent magnetization state m = M^/Mg = cos(6') 
pointing along the z-axis, is not so easily found. One may get this effective tensor as a 
statistical average of the Polder tensor in a non-interacting magnetic grains picture. 
This strategy was investigated^^. Because the distribution of the magnetization 
direction (the magnetic texture) has to be given, and can change as a function of m, it 
results that the non-interacting assumption is strictly satished for a system consisted 
in an assembly of single-crystal spheres situated in a non-magnetic medium far from 
one another, a condition rarely observed in dense soft ferrites2^>^. No interaction 
between domains and domain-walls can be considered in such picture and this also 
neglects the influence of the shape of the magnetic domains and pores. So a two 
steps mechanism is introduced. The magnetostatic interaction between domains has 
to be treated first and then the resulting composite medium has to be statistically 
averaged. 

For any value of m, consider that the effective 3x3 permeability tensor pg can 
be diagonalized and expressed in a circular frame by the following diagonal tensor 
(see Appendix 1X1 for the notations): 




/Xg “1“ Kf, 0 

0 ^g /xg 

^ 0 0 


0 ^ 
0 

^^ez j 


( 1 ) 


where /Xg and Kg are the diagonal and off-diagonal components of the permeability 


4 




expressed in a cartesian framei^. is the permeability value which connects both 
the magnetization and the ac held along the z-axis. 


A. Two-phase model 


By using the same argument to the scale of a magnetic ’’grain”, let us consider 
”up” and ’’down” borderless single domain grains which differ only by the direction 
of their magnetization vector and thus by a change of a sign. Their respective 
permeability tensors are 


hi 


^ ^ + K 0 0 ^ 

0 yU — K 0 
y 0 0 j 


and ^2 


^ ^ — K 0 0 ^ 

0 p + K 0 
y 0 0 j 


( 2 ) 


Now, one consider that each borderless domain can be represented with Polder ex¬ 
pressions in frequency. Such analytical expressions of p, k and /i^ are determined 
from the solution of the Landau-Lifchitz equation in the small signal approximation 
such as : 


T] — taQ 

n 

{r]-iany 


( 3 ) 

( 4 ) 

( 5 ) 


with r] = Hk/Mg the reduced anisotropy held assuming an uniaxial symmetry, = 
oj/j^oMs the reduced pulsation with 7 the gyromagnetic ratio. Mg the saturation 
magnetization, /tq the permeability of free space and a is the damping constant. 

If the domain structure (”up” and ’’down”) is distributed such as one cannot dis¬ 
tinguish any direction in the perpendicular plane to the magnetization on a macro¬ 
scopic scale, then the permeability tensor must be diagonal with respect to such a 
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rotating field^. To get that, one may consider domains in the shape of infinite circular 
cylinders^. The demagnetizing tensor for a single cylinder embedded in a effective 
anisotropic medium in the quasi-static approximation is then 


/_ 1 


_ 2/i, 

r= 0 


( 6 ) 


2p,e 

^ 0 0 Oy 


Even when the medium is characterized by an anisotropic tensor, Eq. (IA12D yields a 
diagonal T tensor with Tn — ^22 for both cylindrical and spherical crystallites. In 
our situation, T is invariant through any cylindrical rotation around the 2 ;-axis and 
preserves this structure in the rotating frame. 

Now some EMA procedure has to be applied in order to get the effective quantities 
(pe, i^e, fJ'ez) functions of (p, K, /i^) and m. It is well known that the EMA depends 
on the choice of the reference permeability^^. Among these choices, the strong- 
coupling interaction between domains is assumed and the symmetric self-consistent 
situation is so retained. In our situation, this corresponds to let equation flAlSj) to 
be zero. For each two-phase grain, the local magnetization of the composite medium 
is m = (ui — V 2 )/{vi + V 2 ) where Vi (resp. ^ 2 ) is the volume occupied by the portion 
of the grain where the magnetization is pointing in the 2 ;-direction (resp. opposite 
to). Because the total volume is V1+V2 = V, one can evaluate Vi = V{l + m)/2 and 


V 2 = V{1 — m)/2. Thus Eq. dAlSD reads 

(1 + m)[T-r(^i + (1 -m)[T-r(^2 -^e)]"^(^2 -%) = (7) 


By substituting matrices ([I]), ([2]) and ([S]) into equation ([7]), Pe, are connected 

with /i, K, ^2 and m as 



( 8 ) 
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(9) 

( 10 ) 


A* 




Remarkably when m = dzl, /ie = /x, = k and jj,ez = 1, which implies that the 

Polder tensor is recovered to the saturation limit. Equation (|8]) and ([9]) agree with 


equations (3a) and (3b) of reference^ which state that for any cylindrically symmetric 


but otherwise arbitrary configuration, this EMA procedure does not depend upon the 


details of the domain conhguration:^'^. Equation ([9]) also agrees with several worksi - 


which relies linearly the static magnetization along the z-axis with the value of the 
off-diagonal permeability, at least for small values of m. 

As a consequence, it is worth noting that for an effective two-components medium 
and any diagonal surface tensor P, the demagnetized situation m = 0 always gives 
Ke = 0 and the following implicit equation for : 

^eh'e r22(h'e) l-^ez') A (1 22{,l-^ei ^^ez^^^^e A r22(h'e) f^ez')^f^ ^ ) h" (H) 

depends only on the single component r 22 = Pii- This result is a generalization of 
Eq.dH]) to any geometry and anisotropic embedded medium resulting in a diagonal 
surface tensor. By inserting Eq.(l6]) into Eq. ffTTl) . the equation (0) is recovered. 

B. Three-phase model 

Neglecting any magnetic after-effect or magnetic viscosity^, the permeability 
spectra are described by two types of magnetizing processes: gyration of domains 
and domain-wall motion^^ which are generally analyzed separately^. In order to 
take into account the lower frequency response of the previously introduced biphasic 
medium, the permeability response for a domain-wall has to be consistently added. 
Its spectral response is commonly assumed to be modeled by an isotropic expres¬ 
sion of permeability^^— that follows from a stiff and dampened 180° domain-wall 
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motion2^ : 



( 12 ) 


Here = Hdw/Mg is the reduced domain-wall resonance field, adw is a damping 
factor and Xdw represents its static susceptibility. To allow the spectral representation 
of such permeability to be reduced to two parameters only, its frequency permeability 
may be approximated as a single uniaxial magnetic domain distribution in the case 
of small damping and if Xdw ~ ^/Vdw In order to let the static susceptibility of 
the domain-wall to be higher than the value of a magnetic domain, it is assumed 
that T]dw returns a smaller value of the anisotropy field in the domain-wall than the 
anisotropy field rj characterizing a domain. 

The goal is to treat on the same footing the dynamics of interacting domains and 
domain-wall motion for a given unsaturated value of magnetization m. An EMA 
is thus constructed with a three components system. ”Up” and ’’down” domains 
are considered supplemented by a domain-wall, each component embedded in an 
anisotropic effective medium with a cylindrical geometry. The homogenization pro¬ 
cedure is depicted for a composite medium made up of uniaxial anisotropy particles 
on Fig. (II]). This gives the following implicit equation for /Ig 



(13) 


i=l 


where JIi and 7 I 2 are the permeability tensors of domains and 7 x 3 = Here Vi 

are the volume fraction of each component. As a consequence of the cylindrical 
geometry, the same surface matrix Eq.([ 6 ]) is taken into account for all the embedded 
media. The fraction volumes have to be normalized according to hmm^.±i U 3 = 0 
and one assumes 


2(l + a(l-m2))’'^^ 


1 -|- m 


2(1 + a(l 


1 — m 
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Figure 1. Schematization of the homogenization procedure for a composite medium made 
up of uniaxial anisotropy domains and domain-wall. 


where a describes some proportion of domain-walls in the whole system. When a = 0 
the previous two-phase situation is recovered and the limit case a —)■ oo provides 
= ^^dw■ Because Fas is here zero, the 2 ; 2 ;-component of equation flT3|) gives exactly 
/iej; and reads 

fiz + andwi^ - 

= l + a(l-m^) ' 

One shows that fiez = /^z when a = 0orm = ±las expected. In the demagnetized 
state m = 0, ^^z is the volume addition of both the permeability of the domain 
and domain-wall by reason of the superposition of the z-component of the magnetic 
induction in all media. When Kg = 0, /Xg is the solution of a single equation in flT^ 
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which is written directly in a third-order polynomial form 
3 ^ - a) -f 2 ap 2 - K ^){1 - a)+ 2 a^ldw^^ ,, 

/^e ~r 1 , /^e , /^e f-^dw 

1 + a 1 + a 




^2) = 0. (16) 


The 3-body composite medium appears isotropic but the value of pe depends sur¬ 
prisingly on the off-diagonal k term. Among all the three complex-valued solutions 
of equation (1T6|) . the unique root with the positive imaginary part is followed con¬ 
tinuously in frequency to get the physical picture of a lossy material. 

For m 7 ^ 0, the xx and j/j/-components of equation flT^ expand as 


(1 -I- m) 
(1 -I- m) 


/i - K - /ie + Ke 

fJjQ — K ^ Ke 

fd K jJj^ 

fJj^ fji -\- K — Kq 


+ (i 
+ (i 


Pe T ^ 

Pe + Ai + ^ T 


. /i K Pg 

m) - 

/ig -h /i - K - Kg 


+ (1 



2 (Z (^fddnj /Xg /xg) 

f^e H“ t^dw H“ 

2 ci /Xg Ke) 

/^e “1“ f^dw 


( 17 ) 

and have to be solved simultaneously to get pg and Kg. These equations are identical 
by interchanging k ■H- —k and Kg -f-)- —Kg which guaranties an exhibition of the 
quantities k^, and KKg. A bi-dimensional root-finding numerical algorithm in the 
complex plane using a variant of the Newton procedure is developed to get the couple 
{fie, ^e) for any values of m, a, fi, k and fidw As for the two-phase situation, Eqs. flThll 
and (1T7|) show that if m = ±1 then fie = fi, Ke = n and /ig^ = 1. This is also true 
by construction of the three-phase model. This implies that /Ig goes to the spectral 
Polder tensor to the saturation limit. 


0 

0 


C. Anhysteretic texture 

The presence of a large amount of such composite media not necessarily dis¬ 
tributed at random, has to be assessed by a space dependent fluctuation in the 
magnetization direction to represent its magnetic texture. The statistical average 
of any function U{d,(f)) of the direction angles 9 and 0 is defined as the following 
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( 18 ) 


weighted normalized integraU^: 

^27r pTT 

{U)= dcf) d0sin0/(0»f/(0,0), 

Jo Jo 

where f{6, 0 ) is a normalized texture function that distribute the angles on the unit- 
sphere at equilibrium and is a function of external parameters such as applied dc 
magnetic held, constrains etc. 

When a full demagnetized state is desired at the macroscopic level, then the 
average remanent held (m) should remain zero but on a statistical distribution of 
the local direction of the magnetization. It is particularly true for the uniform angular 
distribution. It is noted that even if Kg 7 ^ 0, one could also realized (ne) = 0 for any 
angular distribution f{9, 0 ) of the composite medium by satisfying 

(ne) = j j 4>) ^ sin^ 0 sin 20 -|- cos 9^ =0, (19) 

which is true when / is uniform. This means that if desired, a locally anisotropic 
composite medium which behaves statically isotropic may be constructed, once their 
composite media are properly distributed. 

For an uniform distribution of composite domains (2- or 3-boby) in the demagne¬ 
tized state, one has (ne) = 0 and 

(he) = j d(f) j d9 sin 9^ [fie + (/^e^ - he) sin^ 0 cos^ 0 ) = ^ , (20) 

il^ez) = [ [ d9 sin 9^{fiez - = (21) 

Jo Jo 47 r 3 

for biphasic cylindrical inclusions uniformly distributed into space and embedded 

in a gyrotropic ehective medium. So the statistical average ehective tensor of a 

demagnetized ferrite is 

{%){m = 0) = ^Tr(^jT, (22) 

which is the result found by Schlomanrif^. 
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For an unsaturated medium when m increases in the ^-direction, the distribu¬ 
tion of the local magnetization angles cannot be uniform and is distributed around 
the 2 ;-axis. Moreover whatever the value of m may take, = 0 for an uniform 
distribution, which is in contradiction to experiments. 

This problem can be circumvent by recalling that the probability that the local 
magnetization occupies an inhnitesimal solid angle in a given direction is therefore 
strongly dependent on the magnetic free energy in this direction. Let 

^27r niT 

(m) = / dcj) d9 sm9f{9, (j)) cos9, (23) 

Jo Jo 

the statistical average of the remanent magnetization. So when all the collective me¬ 
dia are in the paramagnetic state driven by the exterior constant induction held 
pointing in the z-direction, the function / is then given by an equilibrium Boltzmann 
distribution f{9, (p) = X exp(X cos(6'))/47rsinh(X) where X = AgMs^oH which nat¬ 
urally breaks the space isotropy when H ^ 0^. The local magnetic moments of 
the composite medium do not follow nor instantly nor spatially, the direction of 
the external induction by coherent rotations but instead in defocusing their local 
magnetization about the average direction. The free energy is scaled by the Arm¬ 
strong’s parameter Ag which narrows statistically the distribution / when both the 
nature and degree of disorder are known. An estimation of this parameter is given 
by equating linearly (m) to H for small value of X and one has in our case 


A. 


I (h-ez) 1 


(24) 


hoM2 ’ 

a result also found previously^. 

Now the statistical average of /Ig has to be derived. By regrouping all the terms, 
one can link exactly {fie), (^e) and (fiez) with (m) to read 

(m) 


(he) — he + f^e) 

C \{m)) 


(25) 
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(26) 

(27) 


{He) = {m)Ke 

2{m) 

where C{x) = coth(a;) —l/x is the Langevin function and C~^{x) is its inverse function 
such as C~^oC(^x) = CoC~^(^x) = x. These three equations constitute the central re¬ 
sult of this paper. One recognizes that (m) = C{AsfioMsH) is the well known param¬ 
agnetic result. Moreover the external held is modulated by the Armstrong parameter 
and its contribution to rotate the average magnetization depends on the susceptibility 
of the material. In a general situation, both /ig and have to be evaluated by consid¬ 
ering a direct map between the local remanent state m of the composite medium and 
the macroscopic remanent magnetization (m) in a some sort of mean-held approx¬ 
imation, i.e. m = (m). For practical applications, an approximation of the inverse 
Langevin function is taken^^ and one has (m)/£~^((m)) (1 — (m)^)/(3 — (m)^) for 

all values of (m). One observes that Tr((pe)) = Tr(pe) is an invariant quantity what¬ 
ever the value of (m) may take. Moreover when (m) = 0 then {m)/C~^{{m)) = 1/3, 
(fie) = (2/ie + il-^ez) = (^he + /^ez)/3 and (^e) = 0 as expected. This is 

provided by the fact that 4^) = l/dvr. 

In the previous description of the ehective permeability tensor, the magnetization 
ratio (m) is a hxed parameter that states the magnetization along the ^-axis. How¬ 
ever the ehective anisotropy held, which is characterized by the value of rj is assumed 
to be independent of such a state. When the material is probed by increasing an 
external dc held H, a continuous variation of (m) occurs. Neglecting any supplemen¬ 
tary demagnetizing helds coming from the geometry enforced on the domains and 
domain-walls, the ehective gyromagnetic resonance frequency of the magnetic do¬ 
main shifts linearly with the amplitude of H during this process. This experimental 
phenomenon has to be reproduced and to take it into account, the local ehective held 
inside each uniform media has to be carried out from a magnetization law. Several 



il^ez f^e) 
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strategies have been derived to get an anhysteretic or hysteric magnetization law, 
from coherent rotation mechanisn>i^>^ to more elaborated models^^*^. As a hrst step, 
for each branches along the hysteresis loop, a local equilibrium function f{6,(j)) has 
to be evaluated to derive the hysteretic unsaturated spectral permeabilities. Even 
in the anhysteretic situation, one has at least to replace in the calculation of the 
local permeability tensor of domains and domain-wall, the anisotropy field by an 
effective held that rotates statically the magnetic moments from the easy axis of the 
corresponding medium. This is done in our case by adding H/Mg to the anisotropy 
constants, hence expressed as a function of (m) only one has 


T] ^ T] + 
Vdw t 


3((a«>-1)’ 


(28) 

(29) 


By reason of the dependance of {fiez) the anisotropy helds, equation fl^ has to 
be solved self-consistently. This procedure is initiated hrst by considering known the 
case (m) = 0 as a guess of — 0)) secondly, that an inhnitesimal increment 

of (m) does not ahect the previous value of {fiezi^ = 0)) hence found. This procedure 
generates the desired permeability as a function of (m) by consecutive steps such as 
{fiezi^ = 0, (m) + 6{m))) ^ {fiezi^ = 0, (m))). 


III. RESULTS 

Soft NiZn ferrites have been synthesized and the ehective permeability component 
in a coaxial wave guide at the APC7 standard have been acquired^^. These ferrites 
have been demagnetized by proper thermal treatment above the Curie temperature. 
It has been verihed by measuring hysteresis loops along and perpendicular direction 
to the hollow of the cylinder that the average remanent held (m) does not exceed 
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10 This spectral measurement reported on figure ([2]) allows to fit model param¬ 
eters from equation fl2^ . which are collected on table [B The fitted values {ri,a) 



Figure 2. (Color online) Experimental real part (a) and imaginary part (b) of the effective 
permeability (pe) measured in a coaxial wave guide at the APC7 standard for a demagne¬ 
tized sintered NiZn ferrite (points). The black curve comes from the calculated effective 
permeability described in the text using parameters found into table [B 


for domains are in agreement with reference^. The ratio between domain-wall and 
domain frequencies is ranged from (2 — 5)10“^ in MgFe ferrites^^ which is compatible 
with our fitted values. It has been reported a strong dependance on compositioir^, 
porosity^, grain size^*^ and applied stresses^ in the low frequency regime, where 
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a 

V 

a 

'Hdiv 

^dw 

5.0E-1 

6.6E-3 

2.4E-1 

4.0E-4 

9.6E-1 


Table I. Parameters of the described model for demagnetized NiZn ferrites 


domain-wall motion is the supposed dominant mechanism. However, the reported 
dependancies may be interpreted more as a consequence of the EMA that mimic 
the effective permeability tensor than a change of the intrinsic values of the param¬ 
eters of the magnetic inclusions. For example, the static effective permeability of a 
demagnetized medium is given by 


(h'e(O)) 


1 -|- CL^^dw T (P- pciui)(l T \/( (p Pdui) (1 T 4(1 -|- fl)^pPcijj; 

3(1 -(- cb) 


(30) 

where p = 1 -|- l/r/, ^dw = 1 + ^/Vdw are the corresponding static permeabilities. 
It contains an equal contribution of the domain-walls and uniform spin rotation in 
contradiction to the classical interpretation given by Snoek— which links the initial 
permeability to domain-wall susceptibility only. This formula also contains the frac¬ 
tion of the domain-walls a as a continuous parameter. For large value of a, the initial 
permeability is well described by the permeability of the domain-wall in agreement 
with the classical interpretation. On the other side, when a —)■ 0 this is the static 
part of the domain permeability that drives the initial permeability. The classical 
interpretation has already been questioned^ and the assumption of a perfect addi¬ 
tion of the two types of magnetic processes as described in reference^, reveals itself 
puzzling with a varying remanent held. 

The behavior of the effective diagonal and off-diagonal permeability with an in¬ 
creasing value of the remanent magnetization (m) is now investigated. For the same 
parameters in the table HI the spectra are depicted on hgure [3l The magnetic behav¬ 
ior of these spectra is strongly similar to those measured from reference^, once given 
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Frequency (Hz) 


Figure 3. (Color online) Spectral real (straight line) and imaginary (dotted line) part 
of the effective permeability {^e) (sub-figure a), effective off-diagonal susceptibility (^e) 
(sub-figure b) and effective zz permeability {^iez) (sub-figure c) as a function of several 
macroscopic remanent magnetization values (m). 


the magnetization law that connects (m) to the external nniform magnetic field H. 
Becanse the domain-wall permeability tensor does not carry any off-diagonal expres¬ 
sion, the resnlting effective off-diagonal permeability is small at very low freqnency 
and simply dominated by the snsceptibly coming from the domains. This is in agree¬ 
ment with first direct observation of domain rotation in ferrites^ and by accurate 
measurements given by Green et al^. 
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IV. CONCLUSION 


Effective permeability tensor for unsaturated polycrystalline ferrites are derived 
through an EMA and combines domain-wall motion and rotation of domains in a 
single consistent framework. The dispersion of the local magnetization axis that 
encodes the polycrystalline character is taken into account by averaging the free 
energy to restore the magnetic anhysteretic behavior. The initial permeability is 
given as a mixture formula of magnetic inclusions and fraction of domain-wall which 
gives a picture of the low-frequency permeability spectra as a magnetic scattering 
in geometrically arrangement of domains and domain-walls instead of varying mate¬ 
rial properties. It is envisioned that this theory can be extended to treat multiple 
scattering of electromagnetic helds in such cylindrical geometry to include multiply 
peaked spectra observed in several uniform, polycrystalline materials. 
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Appendix A: EMA formalism 

The double overbars indicate a tensor notation of rank 2, such as the i-th an j-th 
component of a tensor A is [A\ij = Aij. A medium of volume V, bounded by a 
surface S is characterized by a permeability tensor /I as a function of frequency and 
whose values are random in position. It is random in a sense that the conhguration 
average of 7i(r), denoted (Jl}, is independent of r. It is also implicitly supposed 
that this conhguration average is equivalent to a volume average, which is the zero 
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wave-vector of the spatial Fourier transform of 7l(r), such as (p) = V~^ fyjl(r)dr. 
The measurable permeability tensor is the constant of proportionality between 
applied magnetic held Hq such as (H) = Hq and the average induction {B) = 
Consider that the discussion is restricted to the quasi-static limit when d/\ 1; 

d is a typical dimension of the scatterer and A is the wavelength of the incident 
oscillating magnetic held for that only the lowest order of scattering wave mode is 
evaluated for high-resistive ferrite^. In this limit, the magnetic held equations can 
be approximated by magneto-static equations only 


V ■ = 0 

V y. H if) = 0 


(Al) 


which are combined up to a given gauge, to hnd a magnetic potential 0(r) inside the 
medium satisfying 

V ■ (7l(r)V0(r)) = 0. (A2) 


Now a magnetic tensor /Iq is considered to depend on the magnetic held in the average 
medium. In the quasi-static limit, the boundary conditions of a uniform applied held 
are imposed and pg does not exhibit any spatial variations. The local permeability 
tensor, thus decomposed as Jlif) = 'J1q + djlif), leads to the following boundary-value 
problem 


V ■ (7IoV0(fO) = -V ■ ((57l(f)V0(fO) in V, 

4>{f) = —Hq ■ r on S. 

With the introduction of the two-points Green’s function g{f,r') dehned by 


(A3) 


V ■ (/roV( 7 (r, r')) = —5{f—r') in V, 
g{r, r') = 0 r' on S, 


(A4) 


the potential 0(r^ admits a formal solution as 

0(r) = —Hq ■ f— / g{f,r')V' ■ {6'Jl{r')'V'(j){r'))dr' (A5) 

Jv 
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and the magnetic field i7(r) = —V0(r), deriving of it, is given by 


= Hq — Jy{5n{r')H{r'))Vg{f,r')dr', 
= ^0 + jyG{r^r'){5']l{r')H{r'))dr’^ 


(A6) 


with 


G'a/3(r,r') 


d‘^g{f, r') 
dr'^drp 


(A7) 


By combining these eqnations, the problem thns rednces to the task of compnting a 
snsceptibility tensor (y) snch as = x(r)i?o and 





G{f, r')x{r')dr' 


'V 


which finally gives 


he = ho + {x{r)). 


(A8) 


(A9) 


If r lies in a medinm labeled i, of volnme n*, eqnation flASp decomposes itself as 
x(r) = 1 + / G{r,r')x{r')dr' + / G{f,r')xir')dr' | , 

V Jvi Jv-Vi ) 

with — ]Iq, and /Ij is the permeability tensor of the medinm i, assnmed 

spatially nniform and known. The last integral is approximated by 


/ G{r,r')x{r')dr' ^ / G{f,r'){x{r'))dr' (AlO) 

Jv-Vi JV-Vi 

once neglecting higher order correlation terms. This last expression closes the integral 
eqnation for x(r). To see this, this last integral is snbstitnted, integrated by parts 
and by imposing the bonndary condition one obtains 

(T - (All) 

where Xi stands for x(r) for r E Vi. Here Tj is a snrface integral, also called depolar¬ 
ization matrix and contains g{f, r') which goes over to the free-space Green’s fnnction. 
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satisfying the differential equation flA4l) and the boundary condition g{f,r') —)■ 0 as 
||f^— r'W —)■ oo. By reason of the translational invariance, it becomes a function of 
f—r' to finally reduces to a constant in space variables. For example, the field inside 
ellipsoids is uniform and x(r) becomes independent of the position. In such case, 
each cartesian component of this matrix is given by 


pa/3 _ _ 


dg{f 


n'pcpr', 


IS' dr a 

where n'^ is the component of the unit normal outward from the surface S 
taking the average of Xi-, one has 


(A12) 

By 


F. = 7^0 + ((1 - < 7 rr)-'>-‘((l - i/JEj-Vp) (A13) 

where 

((T - S0)-') = Jim i..(T - 3Rr.)-' (A14) 

i 

((T- = Jim (A15) 

i 

Now, a two-phase medium of magnetic objects labeled by i = 1, 2, each fully charac¬ 
terized by their respective position and field independent permeability tensors is 
considered. If 7 I 2 , Pg and Fi 2 = F are all hermitian matrices, and if = JIq, one 
shows that equation (IA13D is equivalent to the Clausius-Mossotti-Maxwell-Garnett 
expressioii^ 


(1 - (/ig - Pi)F) ^(/ig-pi) =/(I - (p2-hi)r) ^(h2-/^i) (A16) 


with / is the volume fraction of object 2 . 
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